Physical mechanisms of interface-mediated intervalley coupling in Si 
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The conduction band degeneracy in Si is detrimental to quantum computing based on spin qubits, 
for which a nondegenerate ground orbital state is desirable. This degeneracy is lifted at an interface 
with an insulator, as the spatially abrupt change in the conduction band minimum leads to inter- 
valley scattering. We present a theoretical study of the interface-induced valley splitting in Si that 
provides simple criteria for optimal fabrication parameters to maximize this splitting. Our work 
emphasizes the relevance of different interface-related properties to the valley splitting. 
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Semiconductor nanostructures based on GaAs and 
Si 0, H are approaching the limit where device func- 
tionality relies on degrees of freedom of individual elec- 
trons. Recent progress in material processing allows pre- 
cise controlled doping, band-structure engineering, and 
fabrication of high quality heterojunctions, which in turn 
pave the way for challenging applications such as the de- 
velopment of a scalable solid state quantum computer. 
The past few years witnessed tremendous experimental 
progress in the study of spin qubits at GaAs/AlGaAs 
quantum dots which raises intriguing questions on 
the feasibility of spin qubits in Si quantum dot [3J] or 
donor states [3] at a Si/barrier-material interface. 

A clear advantage of spin qubits in Si over GaAs is the 
long spin coherence times in Si [5[. On the other hand 
bulk Si conduction band edge is six-fold degenerate, a 
complication not present in GaAs. Near a (001) interface 
with a barrier material, this degeneracy is partially lifted, 
with the interface electron ground state remaining doubly 
degenerate. For electron spin qubits, the residual orbital 
degeneracy is an important spin decoherence source [6|]. 
This effect can be overcome if the ground state degener- 
acy is significantly lifted, which occurs close to an inter- 
face that can efficiently scatter carriers between the two 
degenerate valleys that are near opposite ends of the Bril- 
louin zone Measurements of the doublet splitting, or 
valley splitting, present significant variations among dif- 
ferent Si/barrier samples, ranging from to ~ 1 meV Q. 
In this context, a simple physical model that can help 
identify the relevant fabrication-related parameters in or- 
der to maximize the valley splitting is a valuable tool in 
assisting current experimental efforts. 

Theoretical approaches to describe the electronic be- 
havior in the presence of an interface or heterojunc- 
tion range from the effective mass approximation (EMA) 
0, E3] , tight-binding models 11 1 to first-principles enve- 
lope function approach The present study, based on 



the physically motivated EMA [13[ , aims to identify the 
relevance of sample-dependent parameters to the valley 
splitting. Our approach, while simple, is original and per- 
mits the study of the intervalley coupling due to a single 
interface, not employing periodic boundary conditions. 
The full plane wave expansions of the Bloch functions 
at the two conduction band minima obtained from ab 
initio calculations [3] contain relevant physical informa- 
tion about the underlying Si substrate, and are included 
explicitly. We identify physical mechanisms that control 
the coupling strength pointing to convenient choices for 
the barrier material and the interface quality. We also 
highlight advantages and limitations of the EMA, in par- 
ticular the oscillations of the valley coupling with the 
arbitrary choice of the interface position within a single 
monolayer distance [9(, which is an EMA artifact. 

We consider the (001) Si/barrier system, so that the 
6- valley Si bulk degeneracy breaks into a 2- valley ground 
state and a 4- valley excited state quartet ~ 20 to 30 meV 
above the ground doublet [IH . We write the Hamiltonian 
for this problem as 



H = H + U(z) - -z, 



(1) 



where H$ is the unperturbed bulk Si Hamiltonian. Trans- 
lational symmetry is assumed in the xy plane, parallel 
to the interface, reducing the interface potential to a z- 
dependent profile U(z) [16[. An electric field F along the 
z direction pushes the electron towards the interface. It is 
instructive to first consider an abrupt interface between 
Si and the barrier material; we model it here by taking 
U(z) in Eq. ([T]) to be a simple step potential 



U(z) = U stcp {z) = U G(z - zi), 



(2) 



with U$ representing the barrier height. The usual bar- 
rier material in Si devices are Si02, corresponding to 
Uq ~ 3 eV, and Sii^Ge^ alloys, in which case Uq may 
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FIG. 1: (Color online) Step model potential U step (z) (thick 
lines) and the calculated ground state envelope functions ^(z) 
(thin lines) for the same electric field F = 150 kV/cm. The 
solid lines correspond to Uo = 150 meV and dashed lines to 
Uo = 3 eV, in which case the barrier potential is above the 
vertical scale of the figure. The calculated valley splitting 
depends on Uo, ^(O)! 2 and the evanescent tail ^(z > 0). 



be tuned by controlling the alloy composition. A typical 
value, for x ~ 0.2 to 0.3, would be Uo rj 150 meV. 

Within single- valley EMA [ijj], the lowest energy con- 
duction band eigenfunctions for Eq. ([1]) are written as 
^u(r) = $>(z)e f z u„(r) where k p = ±k z are the Bloch 
wave vectors of the conduction band minima (ho ~ 
0.84 x 27r/ao). The env elop e function &(z) satisfies the 
effective mass equation [16[ 
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FIG. 2: (Color online) (a) Absolute value of the intervalley 
coupling as a function of the electric field intensity for Uo = 
150 meV. The absolute values of the relevant terms in Eq. J5]) 
are also shown, (b) Representation of Vvo = Vs + Ve in 
the complex plane for Uo = 150 meV and F = 150 kV/cm. 
Although Ve does not affect |Wo| much, it rotates Vvo (i-e, 
changes its phase), (c) Same as (a) for Uo = 3 eV. 



the periodic functions u± in terms of plane waves [14J 

0±=*(z)e ± ^^c ± (G) e iG ' r , (5) 

G 

where G are reciprocal lattice vectors. The expression 
for Vvo then reads 



-ft 2 d 2 
2m, dz 2 



U{z) 



F 



*{z) = E*{z), (3) V VO = c + (G)c-(G')5(G x -G' x )5(G v -G' y )I(G z ,G' z ), 



G,G' 



where m z is the longitudinal effective mass for Si. The 
ground state is numerically calculated through a finite 
differences method. Figure[T]gives solutions for the above 
values of Uq. Strictly within the EMA assumption that 
the perturbation potential varies slowly in the length 
scale of the lattice parameter ao, the ground state would 
remain doubly degenerate, since one obtains equivalent 
solutions for k M = ±koz. This assumption is clearly not 
valid in the case of the step potential in Eq. ([2|) which 
is discontinuous at the interface position z — zi, thus 
coupling the originally degenerate Bloch valley states 
\4>±(r)). Also, the electric potential has a discontinu- 
ity in the derivative at z = zi due to the different values 
of the dielectric constant e in the two materials. The 
interface position zj is initially taken to be 0. 

Since {|0+), 1^-)} are well separated in energy from 
the excited states, an effective low energy Hamiltonian is 



(6) 

where the orthonormality of the x and y components is 
used, since there is no perturbation potential along these 
directions. The last term is an integral 



I(G Z) G' Z ) = 



-|-OC:- 



|*(z)| 2 e 4. 



u e(z) - -z 

e 



(7) 
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where Eq is the ground state energy obtained directly 
from Eq. ([3]), and the coupling due to the perturba- 
tion lifts the degeneracy leading to the valley splitting 
2 1 Vvo I- The quantity of interest determining the split- 
ting is Vvo , also called valley-orbit coupling, given here 
by Vvo = {4>+\H\4>-) , which is a complex number. Note 
that Hq gives no contribution to this coupling. We write 



with Q = G z - G' z - 2k . Terms with G and/or G' ^ 
contribute to Eq. ^ with values comparable to the G = 
G' = term alone. Integrating the contribution of the 
step potential by parts, Eq. ([6]) is rewritten as 

Vvo = Vs + V E + V F . (8) 

Here, Vs is the contribution from the I(G Z ,G' Z ) of 
the form ^ e l Q z \^{z)\ 2 6(z)dz = i^i|*(0)| 2 , sim- 
ilar to the effect of a Q-dcpcndent S function cou- 
pling potential. The contributions to Ve, of the form 
Jo°° e"^ z d ^^ dz, involve only the evanescent part of 
the envelope function penetrating into the barrier region. 
We find Vf, the electric field potential contribution, to 
be negligible in all cases considered here, meaning that 
the "kink" in the electric field potential at the interface 
does not couple \<fi±), thus following the standard EMA 
assumption. The role of the electric field in controlling 
Vvo is mainly to modify the envelope function \T/(z). 

Fig. [5]Ja) shows the absolute value of the coupling, 
\Vvo\ and of the \V$\ and \Ve\ terms, for Uq = 150 meV, 
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as a function of applied field. In this case \Vyo\ is well 
described by the ((-function contributions alone. If one is 
also interested in the eigenstates, the phase in Vyo be- 
comes relevant and the contribution Vg plays a role, as 
illustrated in Fig. [2(b). Fig. [21(c) gives the results for the 
couplings dependence on the field for Uq = 3 eV, and in 
this case both contributions affect |Vvo|- For the fields 
considered here, \Vyo \ increases linearly with F 0, [n[ . 
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FIG. 3: (Color online) Variation of the absolute value of the 
valley orbit coupling with barrier height for F = 150 kV/cm. 
The step model of Eq. ([2]) and the smeared (L = 1 monolayer) 
profile of Eq. (|10|l are shown as the solid and dashed lines 
respectively. Both present a maximum for Uo ~ 200 meV. 

The expression for V$ indicates that it is proportional 
to the product [7 |*(0)| 2 . As Uq increases, |\&(0)| 2 de- 
creases, so maximizing \V$\ involves particular and not 
obvious conditions. The same is true for Ve, for which 
only the evanescent tail (z > 0) of &(z) contributes. 
As shown in Fig. [TJ when Uq decreases the z-range 
that makes significant contribution to the integration in- 
creases. However, the highly oscillatory phases e z( ^ z in 
the integrand usually suppress this contribution. A com- 
parison of Fig. H(a) and (c) shows that the larger Uo 
actually produces a much larger \Ve\, even though it cor- 
responds to a more rapidly decaying evanescent envelope. 
Maximizing \Vg\ and \Ve\ independently does not neces- 
sarily maximize the absolute value of the complex sum 
Vs + Ve = Vvo- The net result is that \Vyo\ may be 
similar in magnitude even for materials with Uq a factor 
of 20 apart, as shown in Fig. [2J The general behavior of 
I Vvo | with the barrier height Uq is illustrated in Fig. [3] 
(solid line) for F = 150 kV/cm. As Uq increases from 
0, a sharp rise in |VVo| is obtained up to a maximum 
coupling value around Uq « 200 meV, followed by a slow 
decay: Very low or very high barriers tend to suppress 
the valley coupling. Tuning Uq involves changing the al- 
loy composition in the barrier, as in the case of Sii-^Ge^, 
or changing the barrier material. Such fabrication pro- 
cesses are in principle experimentally feasible. 

We now consider the effect of the interface width, 
which is disregarded in the step model. We use for U(z) 
in Eq. (Q] interface models that are similar to previously 
measured and calculated profiles [l7|). namely 
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FIG. 4: (Color online) (a) Normalized potential profiles for 
the finite-width interface models[Eqs. © and (|10p ]. The step 
potential is also given. For U sa , uss (solid curve) and U exp 
(dashed curve) the same width L = 1 monolayer is taken, (b) 
Variation with the interface width L of the normalized valley 
coupling for F=150kV/cm and Uo = 150 meV . The step po- 
tential corresponds to L = in both models. Reported sample 
measurements are consistent with L w 0.5 — 1.5 monolayers. 
Note that the slightly steeper (7 gauss significantly increases 
| Wo | ■ The shape of these curves is not strongly dependent 
on the value of Uo- 

U cxp (z) = ^ [tanh(z/L) + 1] , (9) 

u sauss (z) = Y Erfc(_z/jL) ' (10) 

where Erfc(a;) is the complimentary error function. Both 
are characterized by a width L, and reproduce U s t e p(z) 
for L = 0. They differ in the asymptotic behavior (re- 
spectively exponential and gaussian). The curves are 
very similar [see Fig, 4(a) for L = 1 monolayer], as quan- 
tified by the RMS deviation with respect to the step po- 
tential; the U exp RMS is only 8% larger than that of the 
t^gauss profile. Yet, values of |VVo| differ by a factor of 
two. This factor may be as large as 3 (for Lsal.3) with 
the same RMS ratio between the profiles. Fig. |3Jb) gives 
the calculated coupling versus L. 

The interface width suppresses the intervalley cou- 
pling, as illustrated by the dashed lines in Fig. [3J The 
sensitivity of the coupling to the functional form of the 
interface potential for the same value of L indicates that 
in real samples the coupling is strongly dependent on the 
type of interface disorder and roughness. Both effects 
contribute to the experimentally observed variation in 
the intervalley scattering among different samples. 

A fundamental limitation of the step potential is the 
sensitivity of the results to the interface position zj 
(taken to be in the results presented so far) on a sub- 
monolayer length scale. This artifact, already obtained 
by Sham and Nakayama , is due to the periodic parts 
of the Bloch functions, which carry information about 
the underlying Si lattice and lead to different interference 
patterns due to the interface perturbation according to 
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its location. As expected, \Vvo\ shows oscillatory be- 
havior with a 1 monolayer period, as seen in Fig. [5l but 
different models lead to different relative phases. The 
more physical model with a finite interface width par- 
tially overcomes this artifact, as the amplitude of the os- 
cillations is reduced for finite L. This is illustrated by the 
squares and circles in Fig. [5l and leads us to believe that 
a realistic model of the interface would largely remove 
this sensitivity to the exact interface location. 



0.24 
0.16 
0.08 
0.00 

0.24 
0.16 
0.08 
0.00 
-0 



(a) U=150meV 



A " A A A 

a □ D ° □ 0°°° n° 
□ □ □□ a n 

°ooo°°o 00 oooooo°° 



A A 

^ o o o o ° 



(b) U =3 eV 



"nop 

>-ooooodoooooooooooobooooo< 



0.00 
z (nm) 



FIG. 5: (Color online) Dependence of the intervalley coupling 
on the positioning of the interface profile relative to the lat- 
tice atomic positions for Uo = (a) 150 meV and (b)3 eV. The 
triangles represent U s t ep [Eq.(j2j], squares U gaU ss [Eq. (Tl0]| ] and 
with L —1 monolayer ~ 0.13 nm. 
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In summary, we have studied the conduction electron 
valley splitting in Si at a single Si/barrier interface us- 
ing EMA. For a sharp interface, the simplicity of EMA 
allows us to identify two main contributions to the valley 
splitting: the electron wave function at the interface and 
the wave function gradient for the evanescent wave in 
the barrier material. We show that an external field (to 
increase the wave function at the interface) and an appro- 
priate barrier potential height (through proper choice of 
barrier material) can help maximize the valley splitting. 

A shortcoming of the sharp interface model is that its 
results are highly sensitive to the interface position. Tak- 
ing into account the finite interface width reduces the sen- 
sitivity to the interface location while also reducing the 
intervalley scattering, because it blunts the sharp inter- 
valley interference. This is illustrated by our result that 
steeper interfaces always favor larger intervalley splitting, 
while smoother profiles tend to reduce it, and may even 
lead to negligible coupling, as demonstrated by the U exp 
profile in Fig. [5jb) . 

We do not expect our EMA results to be quantitatively 
accurate. Many effects are not included, such as strain, 
interface misorientation (lo| . atomic scale disorder, lat- 
eral confinement, or many-body corrections to the valley 
splitting. Nonetheless, the splittings 2| Vyo | on the order 
of 0.5 meV we obtained are in fair agreement with avail- 
able measurements in Si/Si02 and Si/SiGe interfaces Q. 



In conclusion, we have calculated electron valley split- 
ting in Si at a Si/barrier interface. We show that a size- 
able single-particle splitting can be generated by applying 
a proper external field, choosing an optimal barrier ma- 
terial of a suitable potential height and producing sharp 
interfaces. Lateral confinement and many-body correc- 
tions can potentially further increase this splitting. 
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